Towards single-cycle squeezing in chirped quasi-phase-matched optical parametric 

down-conversion 
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We propose a method for generation of a single-cycle squeezed light by parametric down-conversion 
in a chirped quasi-phase-matched nonlinear crystal. We find an exact quantum solution for this 
process valid for an arbitrary parametric gain, and discover an ultrabroadband squeezing in the 
down-converted light with flat squeezing spectrum comprising the full optical octave. We describe 
a scheme for observation of this kind of squeezing using second- harmonic generation as an ultrafast 
correlator. 



PACS numbers: 42.50.Dv, 42.50.Ar, 42. 65. Re 



Squeezed light is one of the central phenomena in 
the modern quantum optics, being on the one hand a 
meso- or even macroscopic object with essentially quan- 
tum properties, and on the other hand a valuable re- 
source for metrology and for quantum information pro- 
cessing. Both the degree of squeezing and the squeez- 
ing bandwidth are important for potential applications. 
At present, the values of squeezing as high as 11,5 dB 
in a band of 100 MHz ^ and 0.3 dB in a band of 2 
GHz [2] have been observed in CW regime. Current ex- 
periments withpulsed squeezed light reach a bandwidth 
in the THz and even tens of THz ^ range. In this 
Letter, we describe a technique for generating a squeezed 
light with ultimate possible squeezing bandwidth com- 
prising all optical spectrum, 250 THz in a realistic exam- 
ple, the degree of squeezing being almost constant within 
this huge bandwidth. 

The proposed technique is based on the process of 
parametric down conversion (PDC) of light in a quasi- 
phase-matched (QPM) periodically poled nonlinear crys- 
tal with linear chirp of the poling frequency, known to be 
able to amplify or generate ultrabroadband optical fields. 
Linearly chirped QPM crystals were shown to produce 
single-cycle biphotons in the low-gain regime of PDC 
[Toj and to provide amplification of classical pul ses with 
high constant gain over a broad bandwidth |ll| - [l3| . In 
particular, in Ref. [l^ a parametric amplification in a 
mid-IR of pulses as short as 75 fs with parametric gain 
over 40 dB has been demonstrated. Here we present a 
quantum theory of PDC in QPM media in the high-gain 
regime and demonstrate that this process can be used for 
generating utrabroadband squeezed light with squeezing 
spectrum comprising the full optical octave. 

We consider the process of collinear non-degenerate 
type-I PDC, where one photon of the pump wave with 
frequency ujp is annihilated to create one photon of the 
signal wave with frequency wq + ^ and one photon of the 
idler wave with frequency iOQ — fl, where ujq ~ i-Op/2. The 
widths of the signal and idler spectra in PDC arc limited 
by the phase-matching condition. In a QPM crystal with 
linear chirp this condition can be satisfied for a broad 
band of frequencies at different spatial positions @. 



We describe the down-converted field in the Fourier 
domain by photon annihilation operators a(51, z), cor- 
responding to annihilation of a photon with frequency 
ujQ + D, at point z. The slowly varying operators b(fl, z) 
are defined by the relations a{fl, z) — h{VL, z) exp[ik{n)z], 
where k{n) is the wave vector at frequency luq + n. For a 
fixed value of $7 > we have two independent quantum 
operators, z) and b^{—il, z). In a QPM crystal with 
the spatial frequency K = Kq — Qz, where C, is the chirp 
parameter, these two operators are coupled via Q 



dz 



(1) 



where k, is the nonlinear coupling coefficient and A(r2) = 
kp~ [fc(il) + /c(— il) + ifo] is the phase mismatch, kp being 
the wave vector of the plane monochromatic pump wave. 
Note, that the phase mismatch is an even function of the 
detuning i7, as usual for a type-I PDC, which implies an 
important symmetry of Eq. ([T|) and its solution. 

To find the solution of Eq. ([1]) we rewrite it in a more 
suitable form, introducing the operators h{U,,z) by 

a{n, z) ^ b{n, ;^)e*MO).g-(A(n).+c^V2+^+V2)^ (2) 

where ip = arg(K) is the phase of the pump field. Next, 
we introduce a new variable x = A(fl)/ ^/^ + z^/^. In the 
variables {il, x) Eq. Q becomes 



db{n,x) 
dx 



— a;6(ri, x) 



(jb\^n,x), 



(3) 



where cr = |K|/-y/C is the new nonlinear coupling parame- 
ter. The system of two linear first-order differential equa- 
tions, Eq. ^ and its Hermitian conjugate with sign in- 
version for f2, is equivalent to one second-order equation 



d%Vt,x) 
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5(r!,x)=0, (4) 



havin g so lutions in the class of parabolic cylinder func- 
tions 14!]. Let us denote two linearly independent solu- 
tions of Eq. Q as (j)i{x) and (/)2(a;) with constant Wron- 
skian W . These two functions can be chosen among var- 
ious pairs of special functions of the mentioned class. We 
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introduce the "reciprocal" functions cj)i{x), i = 1,2, by 
the relation 



(5) 



By construction, the pairs {(j)i{x),(j)i{x)) are solutions of 
the system created by Eq. (|3]) and its Hermitian conju- 
gate with sign- inverted fl. The general solution of this 
system with boundary conditions at x = xq is written as 



b{n,x) 



b\~n,x) 



a 
W 



a 
W 



4>i{x) 02 (a;) 

4>l{xo) (I>2{xq) 
(f)l{x) 4>2{x) 

(t>i{xo) 02(2:0) 



6(f2,xo) 
6t(-l],a:o)|, (6) 

6t(-f],xo) 
6(il,xo) 



01 (x) 4>2{x) 
(t>i{xa) (t)2{xo) 

01(2:0) 02(a;o) 



Equations ^ represent the Bogoliubov transformation of 
the field operators, which is known to produce multimode 
squeezed states of light Is] . 

For practical calculations we choose the functions 
01 (x) and 02 (x) from the family of Whittaker parabolic 
cylinder functions TiJ]: 



0i(x) = A.(xe^"/'), 
02(x) = i?_i_„(-xe-'^/4)^ 



with the corresponding reciprocal functions 

0l(x) = J.l/2g.3V4^^^_^(^g»V4)^ 

02(.t) = iy-^/\r"/^D^,,{~xe-'^/^), 



(7) 



(8) 



where i/ = <t^. The Wronskian of 0i(a;) and 02 (x) is 

Using these functions as the basis, we find from 
Eqs. ^ the transformation of the field operators 6(r2, z) 
from the crystal input at z = to its output at z = L: 

b{n,L) = A{n)b{n,o) + B{n)b''{~n,o), (9) 
P{-n,L) = A{n)b'<{-n,o) + B{n)b{n,o), 



where 



A{n) 



B{n) 



ae 



D,,ixLe'^/^)D^,,i-xoe-'^/^) 

VvD^^^U-XLe-'^'^)D,,~i{xoe''''^) 
^"/^ I?-i-,,(-a;Le-"/*)A.(xoe*"/4) 



(10) 



-A.(xLe*"/4)i^_i_,,(-a;oe-"/") 



-■Ku/2 



while A{Vl) and _B(il) are obtained from A(57) and B{Vl) 
by mutual exchange of 0i (x) and 02 (x) with their recip- 
rocal functions. 



Frequency enters into Eqs. (ITOl) only via the depen- 
dence of XL = A(rj)/VC-hLVC and xq = A(r2)/v^ of 
the phase mismatch, meaning that A{Vl), B{Vl), A{H) 
and BiVl) are even functions of the detuning VL. It will 
be shown elsewhere [3 that A{Vi) = A*{n), B{n) = 
B*{D,), though in general 0„(x) ^ 0^(x). 

In Eqs. (fTO)) we recognize the Green's functions for 
the signal and the idler waves in a chirped QPM crys- 
tal, obtained by a purely classical treatment of the prob- 
lem [l2| . In our approach Eqs. ([SJ describe a unitary 
transformation of the field operators, realized in the non- 
linear medium. The unitarity of the transformation re- 
quires that \A{n)\^ - |B(17)P = 1, and A{n)/B{n) = 
A{—il,)/B{—n). The first relation can be proven [l^ 
with the help of the recurrence relations for the parabolic 
cylinder functions [l3], while the second relation follows 
from the evenness of the functions A{il) and _B(f2). 

Using Eqs. ^ and ^ we obtain the transformation 
of the field operator a(il, z) in the form of one equation: 



a{n,L) = u{n)a{n,o) + v{n)a\-n,o), (ii) 



where 



u{n) 
v{n) 



^^f2)gifc(n)Lgi(A(n)L+cLV2)^ 



(12) 



= iB{n)e 



ifc(n)Lgi(A(0)L+CLV2)+jv 



Note that these two functions of fl are not even any- 
more, because, due to dispersion, the signal and the idler 
waves acquire different values of phase k{il)L when pass- 
ing through the crystal. The properties of the functions 
U{^) and V{^), required by the unitarity of Eq. pT|) . 
follow directly from that of functions A{n) and B{n). 

With the help of asymptotic properties of the 
parabolic cylinder functions, various limiting forms of 
Eqs. (fTU]) can be explored. In the case of low gain, 
v <C 1, we can take the limit v —4- 0, and us- 
ing the relations |l3l -Do(x) = exp[— x'^/4], £)_i(x) = 
'7r/2exp[x^ /4] (l-ierfi[-ix/\/2]), we obtain 




(13) 
(14) 



|erfi 


(1 -1- i)xo 


— erfi 


'{l + i)xL 


} 
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which, together with the phase factor from Eqs. (fT2|) give 
the transformation, obtained in the low-gain regime as a 
perturbative solution of the initial Eq. ([T]) Q. 

Using the solution given by Eq. pTj) . it is straight- 
forward to find the optical spectrum S{uj) of the PDC 
light at the output of the crystal, defined by the relation 
(at(rj, L)a{n', L)) = S{uJo + n)5{n - O'). Employing the 
canonical commutation relations, [6(ri, 0), 6^(1]', 0)] — 
^(5(ri — f2'), we obtain 



(15) 
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This spectrum is shown in Fig. [l]for the 2 cm crystal of 
LiNbOs, pumped at 0.42 /im and quasi-phase-matched to 
produce down converted hght from 0.46 to 0.75 /zm, as in 
Ref. [s^]. The refractive index is obtained from the Sell- 
meier equation for extraordinary wave in LiNbOa. The 
difference from Ref. [8] is in the pump intensity which is 4 
orders of magnitude higher, corresponding to = 0.146. 
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FIG. 1: Optical spectrum of the ultrabroadband squeezed 
light for the 2 cm crystal of LiNbOa, pumped at 0.42 fim. 
The spectrum is symmetric with respect to the frequency ujo, 
as implied by the evenness of 

As well known jl5|, the transformation given by 
Eq. (jlip generates broadband quadrature squeezing in 
the PDC light. For each pair of modes with opposite 
detunings we construct two quadrature operators as [l5l | 

Xiin,L) ^ b{n,L)e''f''^"'^^ +b\-n,L)e-''f''^"'^\ (16) 

where the angle of squeezing ipi^i ^ 

^aTg[U{Q)V{—fl)] determines the orientation of the 

squeezing ellipse. In terms of these quadratures the 
transformation Eq. (jll[) can be rewritten as 



X^{n,L) = e 



exp[±r(r!)]X^(l],0), (17) 



where the upper (lower) sign corresponds to /i = 1 
in = 2), exp[±r(0)] = |J7(0)| ± \V{n)\, and k^{n) = 
[k{n) - k{-n)]/2. It follows from Eq. (H?]) that the 
quadrature X2{fl,L) is squeezed below the standard 
quantum limit, while the conjugate quadrature Xi{fl, L) 
is stretched above that limit. 

The spectra of the quadratures components are defined 
as follows: {Xf,{n, L)Xf,{n' , L)) = Sf,{n)6{n + n'), the 
spectrum of the squeezed quadrature X2 being known as 
the spectrum of squeezing S2{^) = exp[— 2r(f7)]. This 
spectrum is plotted in Fig. [2] for the same crystal and 
experimental settings as in Fig. [T] 

As follows from Fig. [2l an octave-broad squeezing can 
be generated in realistic QPM crystals, widely available 
today. The degree of squeezing can be estimated by con- 
sidering the limit of a; 3> for the parabolic cylinder 
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FIG. 2: Spectrum of squeezing in the high-gain regime for the 
same settings as in Fig. [1] u — 0.146. 



functions, which results in approximation of the modu- 
lus oiU{fl) within the squeezing band by the Rosenbluth 
gain factor [l^, [13: |t^(^)l ~ e'^"; and consequently, 
S2{^) w (e""" - ^/e'^'^" - 1)2. In the low-gain regime 
52 (ri) « 1 at all frequencies. 

The angle of squeezing ij}{Vt,L) is calculated numeri- 
cally and shown in Fig. |3l It spans 6000 rad through 
a bandwidth of 1.5 x lO^^Hz, or on average 0.004 
rad/GHz. This angle is compared to the compensation 
angle 6'o(f^) = A(f7)V(4C) - [k{n) + k[-n)]L/2 found 
for the low-gain regime @] . We see these two angles have 
similar behavior. Therefore, the angle of squeezing can 
be compensated by the same technique as in Refs. @, Q. 




£LI(jJq 

FIG. 3: Angle of squeezing tp{Q,L) in the high-gain regime 
(blue solid line) and the low-gain compensation angle 9o{Q) 
(red circles and dotted line). The crystal is the same as in 
Figs. [T] and [21 The phase of the pump field is taken as zero. 

In the time domain for coherent field the relation 
{SX^{t)SXf,{t')) = 6{t - t'), where SX = X - (X), im- 
poses the standard quantum limit of the quadrature mea- 
surements. In the limiting case of a flat squeezing spec- 
trum equal to e"^'' from to ujq and a constant 'ijj{ft, L), 
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taking a Fourier transform of Eq. pHl) we obtain 

{dX2{t)dX2{t')) = e~^'-~S{t-t'), 



(18) 



where S{t) = — sinc(cjoi) is a delta-like function with 
the width of the order of an optical period at the car- 
rier frequency ujq. Equations (jl8p describe an ultimate 
limit of squeezing in temporal domain with the char- 
acteristic time as short as one optical period, the phe- 
nomenon which we call "single-cycle squeezing" . This 
kind of squeezing would have numerous potential appli- 
cations, for example, in the quantum metrology of ultra- 
short optical processes. The squeezing parameter r will 
be determined by the available parametric gain, which 
can be very high for nanosecond [llj or femtosecond 
optical parametric amplification. 

The ultrabroadband nature of the generated squeezed 
light can be observed using second-harmonic generation 
as an ultrafast correlatorj_similar to the experiments in 
the low-gain regime [18l - l20l |. When the PDC light is 
directed to a thin crystal allowing the second-harmonic 
generation, the field of second harmonic can be written 
as 02 (t) = a2o(t) + €0^(1), where a2o{t) is the vacuum 
field of the second harmonic at the input, a{t) is the 
field of the PDC light, and e is a small parameter, de- 
termining the efficiency of the process. The spectrum of 
the second harmonic S2{(^), determined by the relation 
(02(1^)02 (w')) = Ssh{^)S{uj — Lu'), reads as: 
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FIG. 4: Normalized photon flux "3?(r) (a,c) and normalized 
amplitude quadrature X{t) (b,d) for the coherent component 
of the second harmonic as functions of the delay time for the 
cases of low chirp (a,b), = 1.14 x 10^m~^, and high chirp 
(c,d), C = 5.64 X lO'^m-^ 



where Qs is the central frequency of the signal, while Uq 
and Vq are the values of |?7(f^)| and within the 

squeezing band. 

The measurement of the dependence of the coherent 
component photon flux $(r) = \Ecoh{T)\'^ on the small 
delay r (Fig. S]) would be, thus, a detection of the ul- 
trabroad bandwidth of the squeezed light. Homodyne 
detection of the second harmonic field with the PDC 
pump as a local oscillator [3 will allow one to observe 



g f^o pump as a local oscillator [8|J will allow one to observe 

— / Uc{n)Vc{-n)dn (19) ^-^^ amplitude quadrature X{t) = 2Re{e-'"^ Ecohir)} , 



2 / \v{n)v{uj - Ljp - n)\- dn. 



where Udn) = U{n)e-'^'^^^\ Vdn) = V{n)e-^'^'^^^ are 
the transformation coefficients including the effect of 
compensation by a phase shift 9{^l). The spectrum 
Eq.([T51) contains two components: the coherent one, 
sharply peaked at the pump frequency, and the incoher- 
ent one, spectrally very broad. The field of the coherent 
component in the photon-fiux units reads as 



(20) 

If the compensation angle matches exactly the squeez- 
ing angle 9{il) — ^{VL,L), the total photon flux of the 
coherent component of the second harmonic reaches its 
maximum value. If, in addition to compensation, a 
small delay r is introduced into the signal wave, 9{Vl) — 
ip{fl,L) + rfl, fl > 0, then the coherent component dis- 
appears as soon as r exceeds the inverse of the squeezing 
spectrum width 27r/Ar2, which for the crystal discussed 
above is 4.1 fs. Approximating the squeezing spectrum 
by a rectangle within the squeezing band, we obtain: 



Ecohir) 



[/oVbAf7e~*''=^sinc [AnT/2] , (21) 



oscillating with delay at the optical frequency with an 
envelope determined by the correlation time, and thus 
providing a direct comparison of the correlation time to 
the optical cycle. In Fig. 0] we show the photon flux and 
the quadrature for a low-chirp crystal with a spectrum 
5 times narrower than that considered above, where the 
correlation time exceeds significantly the optical period, 
and the same quantities for the high-chirp crystal, where 
these times are of the same order. 

Ultrabroaband multimode structure of single-cycle 
squeezed light can be also investigated using the methods 
developed for studying multimode quantum frequency 
combs G], squeezed quantum pulses [2lj, and quantum 
fiuctuations in optical solitons |23l . We shall address this 
topic elsewhere [16]. 

In conclusion, we have presented a method for genera- 
tion of a single-cycle squeezed light as an ultimate limit 
of ultrabroadband squeezed light with squeezing compris- 
ing the full optical spectrum. Our scheme is based on 
parametric down-conversion in a linearly chirped quasi- 
phase-matched nonlinear crystal. Our proposal is within 
the reach of current technologies related to femtosecond 
optical parametric amplification in such media. 

We are grateful to Steve Harris for valuable discussions. 
This work was supported by region Nord-Pas-de-Calais 
(France) via project "Campus Intelligence Ambiante" . 
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